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O I Abstract 

Statistical mechanics is generalized on the basis of an additive in- 

ff^ I formation theory for incomplete probability distributions. The incom- 

K*" ' plete normalization J2T=i P^' = 1 is used to obtain generalized entropy 

,^ ■ S = —kJ2T=iP'i^^Pi- 'T^^ concomitant incomplete statistical mechan- 

cn . ics is applied to some physical systems in order to show the effect of 

Q!^ I the incompleteness of information. It is shown that this extensive gen- 

f— ^ ■ eralized statistics can be useful for the correlated electron systems in 

^^ , weak coupling regime. 
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C : 1 Introduction 

o 



If the coin of the Chevalier de Mere^ comes down heads for 49%, tails must 
be 51% because it has only two sides. Since Pascal, this hypothesis has never 
been questioned. 

^In 1654, the Chevalier de Mere, a gamester, proposed to Blaise Pascal the following 
problem : two players of equal skill want to leave the table before finishing their game. 
Their scores and the number of points which constitute the game being given, it is desired 
to find in what proportion they should divide the stakes. In his answer to this question, 
Pascal laid down the principles of the theory of probabilities. 



Let us define an ensemble ^2 of N elements (e.g. number of tosses of 
a coin). Every element has v possible states (sides of coin). A random 
variable (RV) of this ensemble is denoted by ^ of which the value is Xi at 
state i with probability pi. All observed values of ^ constitute an ensemble 
X = {xi,X2, .■■,Xv} with a probability distribution P = {pi,p2, ..,Pv}- If v 
is the number of all the possible values of ^, then x is a complete ensemble 
and ^ is called a complete random variables (CRV)[T|. In the toss of a coin, 
the coin position is a CRV having two values {v = 2) : heads or tails. ^ is 
referred to as independent CRV if all its values are independent (e.g. the 
result of a toss has not any influence on the result of the next toss) and 
exclusive (e.g., heads can not be tails). In this case, P is called a complete 
probability distribution (CPD) for which we have the following postulate : 

The corresponding mathematical framework for calculating P is sometimes 
called Kolmogorov algebra of probability distribution [1, . Eq. ([T)) is the foun- 
dation of all probabilistic sciences using CPD. For the average value of certain 
quantity O of the ensemble Q, we should write 

<0>=j2PiO^. (2) 

1=1 

where Oj is the value of O in the state i. 

It is noteworthy that the sum in these equations is over all possible states. 
Therefore, logically, all statistical theories constructed within Kolmogorov al- 
gebra of CPD should be applied to systems of which all the possible states 
are well-known so that we can count them to carry out the calculation of 
probability or of any quantity. In physics, this requires not only the exact 
hamiltonian and solutions of the equation of motion, but also the mathemat- 
ical tools allowing to treat exactly the known states in phase space. 

However, there are few real systems about which we are capable of treat- 
ing all information. Many systems in physics theory are nothing but some 
isolated and simplified subsystems of the complex and messy world only par- 
tially known. The information needed to specify exactly these subsystems 
is never completely accessible for us. Thus, strictly speaking, the CPD from 
which a completely accessible information can be calculated is no more valid. 



We should address a more general distribution, that is the incomplete prob- 
ability distributions (IPD) due to the incomplete random variables of which 
we do not know all the possible values (i.e. v may be greater or smaller than 
the real number of possible values). CPD is only a special case of IPD|lj. 
Very recently, a possible statistical theory for IPD was proposed |21 El IH E]- 
This statistics based on the incomplete information consideration has been 
applied to a generalization of Boltzmann-Gibbs statistical mechanics (BGS), 
the so called nonextensive statistical mechanics (NSM)[n]. It is shownPJ that 
this nonextensive incomplete statistics (NIS) can indeed overcome some im- 
portant theoretical difficulties encountered in the last CPD version of NSM. 
In addition, the incomplete nonextensive generalization of Fermi-Dirac dis- 
tribution is proved to be the only generalized quantum distribution (among 
several ones) showing the same characteristics as the distribution of strongly 
correlated heavy fermions and so is possible to be applied to this kind of 
systems [4, 5J. 

In this paper, on the other hand, the idea of IPD will be used to generalize 
BGS for the cases where additive information and physical quantities hold. 
We will show that it is possible to use this extensive generalization of BGS 
to describe correlated electrons in the weak-coupling regime. 

2 About incomplete distributions 

A "ideal coin" has only two sides. So we surely have v = 2. But what 
about a coin having thick edge which possibly remains standing without 
being observed? And what is the value of v for a quantum coin? For a 
coin having overlapped heads and tails? For a coin in fractal or chaotic 
phase space? These questions may appear comic. But many real systems 
can surely be regarded as one of these coins, especially the complex systems 
with correlation or effects of correlation that can not be exactly described 
with analytic methods. For this kind of systems, the equation of motion 
becomes incomplete in the sense that some interactions are missing in the 
hamiltonian and that its solution can not yield complete knowledge about 
the system. This incompleteness of knowledge makes it impossible to apply 
Kolmogorov theory because Eq. ((H) and Eq. Q can not be calculated. In this 
case, statistical theories for IPD are necessary. 

The basic assumption of incomplete statisics 2J is that our knowledge 
about physical systems and the relevant probability distributions are incom- 



plete, i.e. YlT=iPi = Q 7^ IP where w is only the number of accessible states 
in phase space. So one can only write J2T=i ^iPi) = 1 where F is certain func- 
tion of Pi. In the case of complete or approximately complete distribution 
(such as in BGS), F is identity function. In my previous work[2j, in order to 
keep the framework of NSM using Tsallis entropy, I proposed F{pi) = pf so 
that 

w 

Ep? = 1, (3) 

where q is Tsallis entropy index [3 El- Since p, < 1, we have to set g G [0, oo]. 
g = should be avoided because it leads to pj = for all states, p^ can be 
called effective probability which allows to relate the parameter q to observed 
results and so to interactions, pi is the 'true' probability which is physically 
useful only when q = 1 for the cases where information is a priori complete. 
According to Eq.Q, Eq.Q should be written as 

w 

<0>=Y.p10.. (4) 

This incomplete normalization is possible whenever the phase space is par- 
tially known or accessible. It is well known that the systems with complex 
interactions show often fractal or chaotic phase space, in which a complete 
calculation of probability becomes in general impossible because the space is 
not integrable or differentiable due to, among others, singular points. In this 
sense, a possible justification of Eq.® can be seen in a work of Tsallis|7^ dis- 
cussing nonadditive energy and probability distributions on fractal supports, 
although at that stage the work was not connected to anomalous normaliza- 
tion like Eq.Q. Considering some simple self-similar fractal structures (e.g. 
Cantor set), one can obtain : 

tiV{oy ^ ^ 

where Vi{k) may be seen as the segments of the fractal structure at a given 
iteration of order k, V{0) a characteristic volume of the fractal structure em- 
bedded in a (i-dimension Euclidean space, df = ^^ is the fractal dimension, 
n the number of segments replacing a segment of the precedent iteration, m 
the scale factor of the iterations and W = n^ the total number of segments 
at the k^^ iteration. If we suppose that the fractal structure with fc — > oo is a 



(i/- dimension phase space containing homogeneously distributed points, the 
complete microcanonical probability distribution at the k^^ iteration should be 
defined as pi = -4t^ = ^^w — where V is the total volume occupied by the 

state points. This distribution obviously sums to one. The problem is that 
V is an indefinite quantity as A; — ;> oo and, strictly speaking, can not be used 
to define exact probability definition. In addition, V is not differentiable (or 
integrable) and contains inaccessible singular points. Thus the probability 
defined above can not be exactly computed. Now if we define pi = -p^ as a 

physical or effective distribution, then we have Z^l^iP/ [Vo/yiO)Y'^^'^ = 1; 
where Vq is a completely accessible and infinitely differentiable support on 
which the calculation of pt is possible. If we choose Vq = V{0), i.e. a d- 
dimension volume containing the (i/-dimension fractal volume, we get Eq.(j3I) 
with q = df/d. The conventional normalization J2i=i Pi = ^ can be recovered 
when df = d. 

The above example is only a case of equiprobable distribution on simple 
fractal structure, but it illustrates very well the possibility that, in complex 
cases, the physical probability may not sum to one and may sum to unity only 
through a kind of power normalization which is in addition consistent with 
the discussions of reference j7j on the mass calculation and the information 
consideration for porous structures. 

As mentioned above, Eq.(jHJ and (J4)) have been applied to NSM established 
with Tsallis entropy 5*5 = —k{J2T=iPi~T^T=iP'i)/i^~l)- As a matter of fact, 
Tsallis entropy implies an information measure in the form ^ _-^ ~ , a non- 
additive generalization of Hartley formula ln(l/pj)P]. In this case, the nor- 
malization oipf is necessary if we want the entropy to have same nonextensive 
properties and variance as generalized Hartley information. Of course, differ- 
ent form of the function F will lead to different statistics. For example, if we 
use Hartley formula as information measure and define expectation value for 
additive entropy and energy with F{pi) satisfying ^^^^^ = a-inx-l-^nZx)^/-/ 
where 7 is an empirical parameter, maximum entropy will lead to the fa- 
mous stretched exponential distribution[8\ pi = ^e~^^^^^'' where Ei is positive 
energy of the system at state i, a and (3 the Lagrange multipliers related re- 
spectively to normalization and energy constraint U = Z^lLi F{pi)Ei. When 
7 = 1, we can recover \ ^^' = ^ as in BGS. 

' ' ox X 

In general, a system becomes nonextensive if it contains interacting parts. 
But in what follows, we will apply incomplete normalization to extensive sys- 



terns. This attempt is partially inspired by the insufficiency of the nonexten- 
sive quantum distributions of NSM to describe correlated electrons in weakly 
coupling regime, as discussed in reference 0- 




Figure 1: g-dependence of the concavity of the extensive generalized entropy 
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3 Extensive incomplete statistics (EIS) 

Now we suppose that the system of interest has A^ interacting elements and 
the information I{N) needed to specify all the elements, as well as the phys- 
ical quantities, are additive (e.g. H = J2j Hj where H is the hamiltonian of 
the compound system and Hj that of the j^^ element). Under these harsh 
conditions, we can postulate [T] : 



1. KV. 



(no information needed if there is only one element) 



2. J(e) = 1 (information unity) 

3. I{N) < I{N + 1) (more information needed with more elements) 

4- -^(riiLi Ni) = YlT=i H^i) (additivity of the information needed to spec- 
ify simultaneously w subsystems containing respectively Ni elements) 

5. I{N) = Iw + YH=iPiI{Ni) (additivity of information measure in two 
steps, where pj = ^ is the probability to find an element in the i*^ sub- 
system and lyj the information needed to determine in which subsystem 
the element will be found. ) 

Only the postulate 5 is different from the conventional form because pi is 
replaced by pi due to IPD. 

As usual, the postulates 1-4 yield Hartley formula I{N) = InN . The 
postulate 5 becomes InN = I^, + Dili pflnNi which can be recast as 



I^ = J2plln{l/p,). (6) 

Now we define an entropy as follows : 

to 

S = kJ2pnn{l/p^) (7) 

which obviously becomes Gibbs-Shannon entropy (Sgs) when g = 1. This 
limit identifies k to Boltzmann constant. It is straightforward to verify that 
S has the same properties as Sgs- 

For microcanonical ensemble, we have 5* = -Inw which decreases with 
increasing q value. In general, AS = S — Sgs < (or > 0) if g > 1 (or 
g < 0) as shown in Fig.(l). This result is consistent with the fact that the 
entropy of a lattice of chaotic maps increases with strong coupling [H]. The 
reader will find later that decreasing q implies increasing coupling for some 
correlated systems. 

For canonical ensemble, the maximum entropy subject to the constraints 
from Eq.Q and (JH) leads to : 

p, = ie-^^' (8) 



with Z = {YlT=i e^'^^^^}^^'^ where Ei is the energy of the system at state i. 
The Lagrange multipher /3 can be determined hj ^ = kp = j; where T is 
the absolute temperature and U the internal energy of the system given by 
U = J2T=iP'i^i- It is easy to verify that U = —-^InZ. 

For grand canonical ensemble, the average number of particles is calcu- 
lated by iV = J2T=iP'i-^i which can be used as a constraint of maximum 
entropy. The partition function is then given hj Z = {YlT=i e~'^^^^^~'^^^^}^^'^ 
which gives 

For quantum particle systems, the usual procedures using Eq. ^, E^ = 
Hkinkjiek and Ni = Y.k{nk)i will lead to 

where fik is the occupation number of the one-particle state k of energy ek 
and n the chemical potential. "-I-" is for fermions and "-" for bosons. 



4 Examples of incompleteness effect 

Now I am presenting examples of the applications of EIS to some simple mod- 
els. The reader should regard them as demonstrations of the incompleteness 
behaviors of physical systems. 

4.1 "Ideal gas" 

The usual calculations give : 

Z={-K[^![!^]3^/2|l/.^ (11) 

U = —NkT (12) 

and 

(a) = AiVfc (13) 

where h is the Planck constant. 



The effect of the incompleteness can be estimated through the energy 
difference 

AU = U - Uo = {- - l)^NkT (14) 

q Z 

which is positive (implying repulsion type interactions) for g < 1 and negative 
(attraction type interactions) for q > 1 {Uq is the internal energy of the 
conventional ideal gas when q = 1). 

The reader may be surprised by the treatment of ideal gas with a the- 
ory for systems including non negligible interactions. From the viewpoint of 
IS, this is just a good example of the philosophy of generalized theory with 
empirical parameter q which is introduced to "absorb" the effects of "unac- 
cessible" interactions or their effects. In this way, interacting systems may 
be mathematically treated as noninteracting or conventional ones including 
only the describable interactions. We are entitled to use, e.g., p^ /2m for the 
total energy of interacting "free particle" , where p is the momentum and m 
the mass of the particle. The extra interaction energy is "absorbed" in the 
parameter q when it has different values from unity. 

4.2 Transport phenomena of ideal gas 

Let W represent the number of particle (of mass m) collisions happening per 
second per unit volume, a usual calculation [TU] will give : 



W = 2n^aJ^ = n'avxl- (15) 

V gvrm ~ 

where n is the particle density and v = J^^ is the most probable speed of a 
particle. Let the mean free path of a particle be denoted by A, the collision 
time r (duration of A) is defined as follows : 

^ V 2W 4na\l kT ' ^ ' 

In this framework, A does not change with respect to BGS case, but r is 
linked to q and increases with increasing q value. This behavior of r can 
affect the electrical conductivity de of metals with free electron model because 

m 4cr V mkT 
which increases with increasing q value. 
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Figure 2: Extensive generalized distribution of blackbody. It can be seen 
that the total emitted energy and the maximal frequency of the emission 
increase when q decreases. 

4.3 blackbody 

The generalized Planck law is given by 



Pq{^,T) 



STihu^ 



oqfShu 



;i8) 



where z/ is the emission frequency and c the light speed. The generalized 
Stephan-Boltzmann law is given by Eq(T) = -%T^ where a is the usual 
Stephan-Boltzmann constant. The distribution Eq. (jl8|) is plotted in Figure 
2 for different q values to show the effect of the incompleteness of information. 
We note that this EIS generalized blackbody is essentially different from the 
nonextensive onejl]. Firstly, the high energy cutoff of nonextensive black- 
body is absent in Eq. p8|) . Secondly, the Einstein emission and absorption 
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coefficients can be easily shown to satisfy -B21 = -B12 in this work as in the 
case of the conventional Bose-Einstein theory. But with nonextensive black- 
body, the ratio B21/B12 varies from zero to infinity according to q value [T^ . 
which possibly opens a door for the interpretation of the behavior of micro 
chaotic laser^^ which still remains a mystery. 

4.4 Ideal fermion gas 

From the EIS fermion distribution given by Eq.()10|). it is easy to verify that 
the Fermi energy e^ at T = in the generalized version is the same as 
in the conventional Fermi-Dirac theory. The zero temperature limit of the 
distribution is therefore not changed. 

If q is not very different from unity and T is not too high, the derivative 
of n does not vanish only when e ~ e^?. In this case, we can use Sommerfeld 
integral ,14J to calculate Fermi energy ep and internal energy for 3D fermion 
gas. The result is : 

..-4[i-g(^)^l, (19) 

and 

where Uq = ^Ne'p is the internal energy of fermion gas at T = 0. The above 
equations show that the decrease of q leads to internal energy increase and 
a drop of Fermi level. This is because that, when g < 1, the particles are 
driven by the repulsion from the lower energy states to higher ones. We can 
see from Eq. ipTH) that. If g = with T > 0, the repulsion would be so strong 
that all states are equally occupied. On the contrary, if g ^ 00, all particles 
are constrained by the attraction to stay at the lowest states, like the case of 
zero temperature distribution. This q effect can be seen in Figure 3 where 
Eq. (fTO|l is plotted for 1-D fermion gas with different q values smaller than 
unity. 

The heat capacity {Cy)q and the magnetic susceptibility Xg of electron 
gas can be given by : 
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Figure 3: Comparison of IS fermion distribution (lines) at various values 
of q with the numerical results (symbols) of Moukouri el al on the basis of 
Kondo lattice t — J model for different coupling constant J[15]. We note 
the flattening of n-drop at Fermi energy ep with decreasing q value. When 
g — !► 0, we have n = 1/2 for all possible states. When q — > cx), we have n = 1 
for all states below ep and n = for all states above ep for any temperature. 
In the calculations, the density of electrons is chosen to give kf^^ = O.SSvr in 
the first Brillouin zone. We see that the lines fit well the numerical results 
for about J < 2. When the coupling is stronger, a long tail in the KLM 
distributions begins to develop at high energy and can not be reproduced by 
IS distribution. A new Fermi surface starts to appear and a sharp n drop 
(cutoff) takes place at about k = O.Tvr. These strong correlation effects are 
absent in the present IS fermion distribution. 
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and 

where 70 is the conventional coefficient of the electronic heat capacity and 
Xo the conventional susceptibility of electron gas at T = 0. 

With the help of Eq. pHl . the ratio of the effective mass rrith to the mass 
m of an electron can be related to the parameter q as follows : 

mth _ {Cy)g{observed) _ 1_ 
m Cv{theoretical) g^ 

For heavy fermions, we surely have g < 1. The reader will find below that 
this fact is conffimed by the q values found for correlated electrons. We 
remember that above relationships are valid only for q and T values which 
keep the sharp n drop at ep. 

On the other hand, in Figure 3, the calculations are valid for any q value. 
For small q, we note a flattening of the 77,-drop at ep with decreasing q values. 
This is just a behavior of correlated electrons observed in some experimental 
and numerical results [TH f IT!) ! \T7 \ IT ^ \W \ 1^ . Figure 3 shows a comparison of 
the momentum distribution (lines) given by Eq. ljlOj) for 1-D fermions with 
numerical results (symbols) given by numerical simulation based on Kondo 
lattice t — J model (KLM)[16j. We see that IS momentum distribution fits 
well the numerical results for about J < 2, the weak coupling regime [T6]. 
When the coupling is stronger (J > 2), a long tail in the KLM distributions 
begins to develop at high energy and can not be fitted with present fermion 
distribution. In addition, in the strong coupling regime, a new Fermi surface 
starts to appear and a sharp n drop (cutoff) takes place at about k = O.Ttt. 
These strong correlation effects are absent in the present EIS fermion distri- 
bution. Similar result is obtained with other numerical results[TK|l21j. This 
is consistent with the fact that high energy cutoff does not exist in EIS. On 
the other hand, the sharp n drop at a higher ep than the conventional one 
is indeed observed in NIS fermion distribution 0. This result con&ms our 
conjecture that EIS is only valid for weak interaction cases due to its ex- 
tensive nature and that nonextensive effect should be considered whenever 
interactions become stronger. 
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4.5 Conclusion 

From above results, we note that the introduction of parameter q is equivalent 
to replacing k in the conventional BGS theory by kq = k/q. One may say 
that the incomplete g-normalization turns out to be a renormalization of the 
constant without any change of physics. It should be noticed that the value 
of Boltzmann constant has been determined for noninteracting systems. So 
any modification of this value undoubtedly implies physical changes in the 
studied system as well as in the corresponding theory. The modification given 
by EIS is a simple one as we guessed considering its validity for the cases 
where nonextensive effect of interactions can be neglected. The successful 
fitting of the heavy fermion distributions shows that this modification refiects 
at least a part of the reality and can be useful. 

Summing up, the conventional BGS statistical mechanics is generalized 
for extensive systems on the basis of the idea that we sometimes can not know 
all the informations about complex systems so that the physical or observable 
probabilities become incomplete and do not sum to one. A generalized addi- 
tive entropy is obtained by using a so called incomplete normalization with 
a empirical parameter q which is intended to "absorb" the effect of complex 
correlations and can be related to the energy of the studied system. So, in 
general, q < 1 and g > 1 imply respectively repulsive and attractive effect 
of the complex correlations. This extensive incomplete statistics is shown to 
be able to reproduce very well the quantum distributions of correlated heavy 
electrons in weak coupling regime. On the other hand, EIS fails to describe 
strongly correlated conduction electrons and localized /-electrons for which 
nonextensive effect should be taken into account. 



5 Acknowledgments 

I acknowledge with great pleasure the useful discussions with Professors A. 
Le Mehaute, Dr. L. Nivanen and Prof. J. P. Badiali on some points of this 
work. 



References 

[1] A. Renyi, Calcul de probabilite, {Paris, Dunod, 1966)P522. 



14 



[2] Q.A. Wang, Chaos, Solitons & Fractals, 12(2001)1431, for erratum, see 
|cond-ma t/0009343| 



[3] Q.A. Wang, Nonextenisve statistics and incomplete information. Euro. 
Phys. J. B, 26(2002)357 

[4] Q.A. Wang, Quantum gas distributions prescribed by factorization hy- 
pothesis of probabihty. Chaos, Solitons & Fractals, 14(2002)765 

[5] Q.A. WANG and A. Le Mehaute, Extensive form of equihbrium nonex- 
tensive statistics, J. Math. Phys, 43(2002)5079-5089 

[6] C Tsalhs, Chaos, Solitons and Fractals, 6(1995)539; 

C. Tsalhs, J. Statis. Phys., 52(1988)479; 

EMF. Curado and C. Tsalhs, J. Phys. A: Math. Cen., 24 (1991)169; 

C. Tsalhs, R.S. Mendes and A.R. Plastino, Physica A, 261(1999)534; 

Silvio R.A. Salinas and C Tsalhs, Brazilian Journal of Physics (special 
issue: Nonextensive Statistical Mechanics and Thermodynamics), 
29(1999) 

[7] C. Tsalhs, Fractals, 3(1995)541 

[8] R Jund, R. Julien and I. Campbeh, Phys. Rev. E, 63(2001)036131 

[9] A.M. Batista and R.L. Viana, Physica A, 308(2002)125-134 

[10] K. Huang, Statistical Physics, (Wiley, New York, 1987), Chap.5, 2nd 
ed. 

[11] Q.A. Wang, and A. Le Mehaute, Phys. Lett. A, 237(1997)28; 

Q.A. Wang, L. Nivanen and A. Le Mehaute, Physica A, 260(1998)490- 
498. 

[12] Q.A. Wang, and A. Le Mehaute, Phys. Lett. A, 242(1998)301-306 

[13] C Gmachl, F. Capasso, E.E. Narimanov, J.U.Nokel, A.D. Stone, J. 
Faist, D.L. Sivco, A.Y.Cho, Science, 280(1998)1493-1656 

[14] Donald A. McQuarrie, Statistical Mechanics, (New York, Harper and 
Row, 1973)P167 

15 



[15] R. Eder, O. Rogojanu, and G.A. Sawatzky, Many-body band structure 
and Fermi surface of the Kondo lattice, cond-niat/971 1248| 



[16] S. Moukuri and L.G. Caron, Phys. Rev. B, 252(1995)R15723; 

S. Moukuri and L.G. Caron, Fermi surface of the one- dimensional Kondo 
Lattice Model, |i cond-m at/9607T23| 

[17] C. Grober and R. Eder, Fermiology of a ID heavy electron model, 
|cond-mat/9812345 



[18] F. Ronning, C. Kim, D.L. Feng, D.S. Marshall, A.G. Loeser, L.L. Miller, 
J.N. Eckstein, I. Bozovic, Z.X. Shen, Photoemission evidence for a Rem- 
nant Fermi surface and d- wave-like dispersion in insulting Ca2Cu02Cl2, 
Science, 282(1998)2067 

[19] W.O. Putikka, M.U. Luchini and R.R.P. Si ngh, Violation of Lu ttinger's 
theorem in the two-dimensional t- J model, jcond- mat/9803140 



[20] W.O. Putikka, M.U. Luchini and R.R.P. Singh, Does the two- 
dimensional t-J model have hole pockets? _cond-mat/9803141_ 

[21] Q.A. WANG, Correlated electrons and generalized statistics, to appear 
in Euro. Phys. J. B, (2003), cont-mat/0205038 



16 



